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The number 7T (/pal/) is a mathematical constant. It is defined as the ratio of a circle's circumference to its 
diameter, and it also has various equivalent definitions. It appears in many formulas in all areas of 
mathematics and physics . It is approximately equal to 3.14159. It has been represented by the Greek letter 
"7T" since the mid-18th century, and is spelled out as "pi". It is also referred to as Archimedes' constant. 

Being an irrational number, 71 cannot be expressed as a common fraction, although fractions such as 22/7 are 
commonly used to approximate it. Equivalently its decimal representation never ends and never setdes into a 
permanently repeating pattern. Its decimal (or other base) digits appear to be randomly distributed, and are 
conjectured to satisfy a specific kind of statistical randomness. It is known that 71 is a transcendental number: 
it is not the root of any polynomial with rational coefficients. The transcendence of 7T implies that it is 
impossible to solve the ancient challenge of squaring the circle with a compass and straightedge. 

Ancient civilizations, including the Egyptians and Babylonians, required fairly accurate approximations to 7T 
for practical computations. Around 250 BC the Greek mathematician Archimedes created an algorithm to 
approximate 7T with arbitrary accuracy. In the 5th century AD, Chinese mathematics approximated 7T to 
seven digits, while Indian mathematics made a five-digit approximation, both using geometrical techniques. 
The first exact formula for 7T, based on infinite series, was discovered a millennium later, when in the 14th 
century the Madhava-Leibniz series was discovered in Indian mathematics.^ 1 ^ The invention of calculus 
soon led to the calculation of hundreds of digits of 7T, enough for all practical scientific computations. 
Nevertheless, in the 20th and 21st centuries, mathematicians and computer scientists have pursued new 
approaches that, when combined with increasing computational power, extended the decimal representation 
of 7T to many trillions of digitsThe primary motivation for these computations is as a test case to 
develop efficient algorithms to calculate numeric series, as well as the quest to break records.®^ The 
extensive calculations involved have also been used to test supercomputers and high-precision multiplication 
algorithms. 

Because its most elementary definition relates to the circle, 7T is found in many formulae in trigonometry and 
geometry , especially those concerning circles, ellipses, and spheres. In more modern mathematical analysis, 
the number is instead defined using the spectral properties of the real number system, as an eigenvalue or a 
period, without any reference to geometry. It appears therefore in areas of mathematics and the sciences 
having little to do with the geometry of circles, such as number theory and statistics, as well as in almost all 
areas of physics . The ubiquity of 71 makes it one of the most widely known mathematical constants both 
inside and outside the scientific community. Several books devoted to 7T have been published, and record¬ 
setting calculations of the digits of 7T often result in news headlines. Adepts have succeeded in memorizing 
the value of 7T to over 70,000 digits. 
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Fundamentals 



Name 


The symbol used by mathematicians to represent the ratio of a circle's circumference to its diameter is the 
lowercase Greek letter 71, sometimes spelled out as pi, and derived from the first letter of the Greek word 
perimetros, meaning circumference.^ In English, 7T is pronounced as "pie" (/pal/ PY).® In mathematical 
use, the lowercase letter 71 is distinguished from its capitalized and enlarged counterpart which denotes a 
product of a sequence, analogous to how £ denotes summation. 

The choice of the symbol 7T is discussed in the section Adoption of the symbol n. 


Definition 


71 is commonly defined as the ratio of a circle's circumference C to 
its diameter d 


C 



The ratio C/d is constant, regardless of the circle's size. For 
example, if a circle has twice the diameter of another circle it will 
also have twice the circumference, preserving the ratio C/d. This 
definition of 7T implicitly makes use of flat (Euclidean) geometry; 
although the notion of a circle can be extended to any curved (non- 
Euclidean) geometry, these new circles will no longer satisfy the 
formula 7T = C/dP^ 


Here, the circumference of a circle is the arc length around the 
perimeter of the circle, a quantity which can be formally defined 
independently of geometry using limits, a concept in calculus. 

For example, one may directly compute the arc length of the top half 

2 2 

of the unit circle, given in Cartesian coordinates by the equation X + y 



The circumference of a circle is 
slightly more than three times as 
long as its diameter. The exact ratio 
is called n. 


- 1, as the integral :t 11 ] 


7T = 



An integral such as this was adopted as the definition of 7T by Karl Weierstrass, who defined it directly as an 
integral in 1841.^ 

Definitions of 7T such as these that rely on concepts of the integral calculus are no longer common in the 
literature. Remmert 2012, Ch. 5 explains that this is because in many modern treatments of calculus, 
differential calculus typically precedes integral calculus in the university curriculum, so it is desirable to 
have a definition of 7T that does not rely on the latter. One such definition, due to Richard Baltzer,^ 12 ^ and 
popularized by Edmund Landau/ 13 ^ is the following: 7T is twice the smallest positive number at which the 
cosine function equals O^ 9 ^ 11 ^ 14 ^ The cosine can be defined independently of geometry as a power 
series,^ or as the solution of a differential equation.^ 

In a similar spirit, 7T can be defined instead using properties of the complex exponential, exp z, of a complex 
variable Z. Like the cosine, the complex exponential can be defined in one of several ways. The set of 
complex numbers at which exp Z is equal to one is then an (imaginary) arithmetic progression of the form: 


{..., —27 ri, 0, 27 ri, 4iri ,...} = {2ivki | k 6 Z} 

























and there is a unique positive real number 7T with this propertyJ 11 ^ 16 ^ A more abstract variation on the same 
idea, making use of sophisticated mathematical concepts of topology and algebra, is the following 
theorem: ^ there is a unique (up to automorphism) continuous isomorphism from the group R/Z of real 
numbers under addition modulo integers (the circle group) onto the multiplicative group of complex 
numbers of absolute value one. The number 71 is then defined as half the magnitude of the derivative of this 
homomorphism^ 18 ^ 

A circle encloses the largest area that can be attained within a given perimeter. Thus the number n is also 
characterized as the best constant in the isoperimetric inequality (times one-fourth). There are many other, 
closely related, ways in which 7T appears as an eigenvalue of some geometrical or physical process; see 
below. 


Irrationality and normality 

22 

71 is an irrational number, meaning that it cannot be written as the ratio of two integers. Fractions such as -y 
355 

and jy are commonly used to approximate 7T, but no common fraction (ratio of whole numbers) can be its 

exact value.Because 71 is irrational, it has an infinite number of digits in its decimal representation, and it 
does not settle into an infinitely repeating pattern of digits. There are several proofs that 7T is irrational; they 
generally require calculus and rely on the reductio ad absurdum technique. The degree to which 7T can be 
approximated by rational numbers (called the irrationality measure) is not precisely known; estimates have 
established that the irrationality measure is larger than the measure of e or In 2 but smaller than the measure 
of Liouville numbersJ 20 ^ 

The digits of 7T have no apparent pattern and have passed tests for statistical randomness, including tests for 
normality; a number of infinite length is called normal when all possible sequences of digits (of any given 
length) appear equally oftenJ 21 ^ The conjecture that 71 is normal has not been proven or disprovenJ 21 ^ 

Since the advent of computers, a large number of digits of 71 have been available on which to perform 
statistical analysis. Yasumasa Kanada has performed detailed statistical analyses on the decimal digits of 71 
and found them consistent with normality; for example, the frequencies of the ten digits 0 to 9 were 
subjected to statistical significance tests, and no evidence of a pattern was foundJ 22 ^ Any random sequence 
of digits contains arbitrarily long subsequences that appear non-random, by the infinite monkey theorem. 
Thus, because the sequence of 7T's digits passes statistical tests for randomness, it contains some sequences 
of digits that may appear non-random, such as a sequence of six consecutive 9s that begins at the 762nd 
decimal place of the decimal representation of TlJ 23 ^ This is also called the "Feynman point" in mathematical 
folklore, after Richard Feynman, although no connection to Feynman is known. 


Transcendence 


In addition to being irrational, more strongly 71 is a transcendental number, which means that it is not the 
solution of any non-constant polynomial equation with rational coefficients, such as 


5 3 

jr_ _ xf , 

120 6 * 


0.[24][b] 


The transcendence of 71 has two important consequences: First, 7T cannot be expressed using any finite 

combination of rational numbers and square roots or n-th roots such as V31 or VlO. Second, since no 
transcendental number can be constructed with compass and straightedge, it is not possible to "square the 
circle". In other words, it is impossible to construct, using compass and straightedge alone, a square whose 
area is exactly equal to the area of a given circle.^ Squaring a circle was one of the important geometry 



problems of the classical antiquity/ 26 ] Amateur mathematicians in 
modern times have sometimes attempted to square the circle and 
sometimes claim success despite the fact that it is mathematically 
impossible/ 27 ] 

Continued fractions 

Like all irrational numbers, 7T cannot be represented as a common 
fraction (also known as a simple or vulgar fraction), by the very 
definition of "irrational number" (that is, "not a rational number"). 
But every irrational number, including 7T, can be represented by an 
infinite series of nested fractions, called a continued fraction: 

1 

7T = 3 H- 



7 + 


Because n is a transcendental 
number, squaring the circle is not 
possible in a finite number of steps 
using the classical tools of compass 
and straightedge. 


15 + 


1 + 


292 + 


1 + 



1 +- 

1 + 

Truncating the continued fraction at any point yields a rational 
approximation for 7T; the first four of these are 3, 22/7, 333/106, and 
355/113. These numbers are among the most well-known and widely 
used historical approximations of the constant. Each approximation 
generated in this way is a best rational approximation; that is, each is 
closer to 7T than any other fraction with the same or a smaller 
denominator/ 28 ] Because 71 is known to be transcendental, it is by 

definition not algebraic and so cannot be a quadratic irrational. Therefore, 7T cannot have a periodic 
continued fraction. Although the simple continued fraction for 7T (shown above) also does not exhibit any 
other obvious pattern/ 29 ] mathematicians have discovered several generalized continued fractions that do, 
such as/ 3 °] 


The constant n is represented in this 
mosaic outside the Mathematics 
Building at the Technical University 
of Berlin. 





















7T 


= 3 + 


1 + 


6 + 


2 + 


6 + 


2 + 


6 + 


2 + 


6 + 


2 + 


6 + 


2 + 


1 + 


3 + 


5 + 


9+ •• 


Approximate value and digits 


Some approximations of pi include: 


■ Integers: 3 

OO QQQ Qrr rni rn 

■ Fractions: Approximate fractions include (in order of increasing accuracy) -y, jQg, jjg, 16604 , 

S 9 , 3 , 3 ° 32 3 i 5 8 . and (List is selected terms from OEIS: A063674 and 

OEIS: A063673.) 

■ Digits: The first 50 decimal digits are 

3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37510..J 31 l (see 
OEIS: A000796) 


Digits in other number systems 

■ The first 48 b inary (base 2) digits (called bits) are 

11.0010 0100 00111111 0110 1010 1000 1000 1000 01011010 0011... (see OEIS: A004601) 

■ The first 20 digits in hexadecimal (base 16) are 3.243F 6A88 85A3 08D3 1319...t 32 l (see 
OEIS: A062964) 

■ The first five sexagesimal (base 60) digits are 3;8,29,44,0,47[ 33 1 (see OEIS: A060707) 


Complex numbers and Euler's identity 

Any complex number, say Z, can be expressed using a pair of real numbers. In the polar coordinate system, 
one number (radius or r) is used to represent z's distance from the origin of the complex plane and the other 
(angle or cp) to represent a counter-clockwise rotation from the positive real line as follows:^ 


z = r • (cos <p + i sin <p ), 



. .2 

where z is the imaginary unit satisfying z = -1. The frequent 

appearance of n in complex analysis can be related to the behaviour 

of the exponential function of a complex variable, described by 

Euler's formula: 

e lip = cos (p + i sin </?, 

where the constant e is the base of the natural logarithm. This 
formula establishes a correspondence between imaginary powers of 
e and points on the unit circle centered at the origin of the complex 
plane. Setting (p = n in Euler's formula results in Euler's identity, 
celebrated by mathematicians because it contains the five most 
important mathematical constants :^ 35 ^ 36 ^ 

e in + 1 = 0. 

There are n different complex numbers Z satisfying Z n - 1, and 
these are called the "rz-th roots of unity".^ 37 ^ They are given by this 
formula: 



The association between imaginary 
powers of the number e and points 
on the unit circle centered at the 
origin in the complex plane given by 
Euler's formula 


2mk/n (A; = 0, 1,2, ... ,71 — 1). 


History 


Antiquity 

The best-known approximations to 7T dating before the Common Era were accurate to two decimal places; 
this was improved upon in Chinese mathematics in particular by the mid-first millennium, to an accuracy of 
seven decimal places. After this, no further progress was made until the late medieval period. 

Based on the measurements of the Great Pyramid of Giza (c. 2560 BC) some Egyptologists have claimed 
that the ancient Egyptians used an approximation of 71 as y- from as early as the Old KingdomJ 38 ^ 39 ^ This 

claim has met with scepticism J 40 ^ 41 ^ 42 ^ 43 ^ 44 ^ The earliest written approximations of 7T are found in 
Babylon and Egypt, both within one per cent of the true value. In Babylon, a clay tablet dated 1900- 
1600 BC has a geometrical statement that, by implication, treats 7T as y = 3.125.^ In Egypt, the Rhind 
Papyrus , dated around 1650 BC but copied from a document dated to 1850 BC, has a formula for the area of 
a circle that treats 71 as (y) 2 ~ 3.16J 45 ^ 

Astronomical calculations in the Shatapatha Brahmana (ca. 4th century BC) use a fractional approximation 
of « 3.139 (an accuracy of 9xlO _4 )J 46 -' Other Indian sources by about 150 BC treat 7T as 
VlO * 3.1622. [47] 

Polygon approximation era 

The first recorded algorithm for rigorously calculating the value of 7T was a geometrical approach using 
polygons, devised around 250 BC by the Greek mathematician Archimedes. ^ This polygonal algorithm 
dominated for over 1,000 years, and as a result n is sometimes referred to as "Archimedes' constant"J 49 ^ 
Archimedes computed upper and lower bounds of 7T by drawing a regular hexagon inside and outside a 
























circle, and successively doubling the number of 
sides until he reached a 96-sided regular 
polygon. By calculating the perimeters of these 

223 22 

polygons, he proved that y- < n < — (that is 

3.1408 < 7T < 3.1429). [50] Archimedes' upper 
22 

bound of y may have led to a widespread 

22 

popular belief that 7T is equal to —! 514 Around 


n can be estimated by computing the perimeters of 
circumscribed and inscribed polygons. 





150 AD, Greek-Roman scientist Ptolemy, in his 
Almagest, gave a value for 7T of 3.1416, which 

he may have obtained from Archimedes or from Apollonius of Perga! 52 ^ 53 ] Mathematicians using 
polygonal algorithms reached 39 digits of 7T in 1630, a record only broken in 1699 when infinite series were 
used to reach 71 digits.'^ 


In ancient China, values for 7T included 3.1547 (around 1 AD), VlO (100 

142 

AD, approximately 3.1623), and (3rd century, approximately 

3.1556)! 55 ] Around 265 AD, the Wei Kingdom mathematician Liu Hui 
created a polygon-based iterative algorithm and used it with a 3,072-sided 
polygon to obtain a value of 71 of 3.1416^ 56 ^ 57 ^ Liu later invented a faster 
method of calculating 71 and obtained a value of 3.14 with a 96-sided 
polygon, by taking advantage of the fact that the differences in area of 
successive polygons form a geometric series with a factor of 4! 56 ] The 
Chinese mathematician Zu Chongzhi, around 480 AD, calculated that 

■aqr: 

3.1415926 < 7T < 3.1415927 and suggested the approximations 71 ~ 

= 3.14159292035... and 7T ~ y = 3.142857142857..., which he termed the 


Milii ("close ratio") and Yuelii ("approximate ratio"), respectively, using Liu 
Hui's algorithm applied to a 12,288-sided polygon. With a correct value for 
its seven first decimal digits, this value of remained the most accurate 
approximation of 7T available for the next 800 years! 58 ] 



Archimedes developed the 
polygonal approach to 
approximating n. 


The Indian astronomer Aryabhata used a value of 3.1416 in his Aryabha tiya (499 AD)! 59 ] Fibonacci in 
c. 1220 computed 3.1418 using a polygonal method, independent of Archimedes! 60 ] Italian author Dante 

/ 9 

apparently employed the value 3+yr ~ 3.14142! 60 ] 


The Persian astronomer Jamshld al-Kash! produced 9 sexagesimal digits, roughly the equivalent of 16 
decimal digits, in 1424 using a polygon with 3 X 2 28 sides! 61 ^ 62 ] which stood as the world record for about 
180 years! 63 ] French mathematician Francois Viete in 1579 achieved 9 digits with a polygon of 3 X 2 17 
sides! 63 ] Flemish mathematician Adriaan van Roomen arrived at 15 decimal places in 1593! 63 ] In 1596, 
Dutch mathematician Ludolph van Ceulen reached 20 digits, a record he later increased to 35 digits (as a 
result, 71 was called the "Ludolphian number" in Germany until the early 20th century)! 64 ] Dutch scientist 
Willebrord Snellius reached 34 digits in 1621! 65 ] and Austrian astronomer Christoph Grienberger arrived at 
38 digits in 1630 using 10 40 sides! 66 ] which remains the most accurate approximation manually achieved 
using polygonal algorithms! 65 ] 


Infinite series 
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The calculation of 71 was revolutionized by the 
development of infinite series techniques in the 16th and 
17th centuries. An infinite series is the sum of the terms 
of an infinite sequence .^ Infinite series allowed 
mathematicians to compute 7T with much greater 
precision than Archimedes and others who used 
geometrical techniques' 67 ^ Although infinite series were 
exploited for 7T most notably by European 
mathematicians such as James Gregory and Gottfried 
Wilhelm Leibniz, the approach was first discovered in 
India sometime between 1400 and 1500 ADJ 68 ^ 69 ^ The 
first written description of an infinite series that could be 
used to compute 71 was laid out in Sanskrit verse by 
Indian astronomer Nilakantha Somayaji in his 
Tantrasamgraha, around 1500 ADJ 70 ^ The series are 

presented without proof, but proofs are presented in a later Indian work, YuktibhaSa, from around 1530 AD. 
Nilakantha attributes the series to an earlier Indian mathematician, Madhava of Sangamagrama, who lived 
c. 1350 - c. 1425.t 7 °] Several infinite series are described, including series for sine, tangent, and cosine, 
which are now referred to as the Madhava series or Gregory-Leibniz seriesJ 70 ^ Madhava used infinite series 
to estimate 71 to 11 digits around 1400, but that value was improved on around 1430 by the Persian 
mathematician Jamshld al-Kashl, using a polygonal algorithm.^ 71 ] 


Comparison of the convergence of several 
historical infinite series for n. S n is the 
approximation after taking n terms. Each 
subsequent subplot magnifies the shaded area 
horizontally by 10 times, (click for detail) 


The first infinite sequence discovered in Europe was an infinite product 
(rather than an infinite sum, which are more typically used in 7T calculations) 
found by French mathematician Frangois Viete in 1593;f 73 ^ 74 ^ 75 ^ 


2 _ y /2 + ][ 2 + V 2 + 

7T “ ~2 2 2 

The second infinite sequence found in Europe, by John Wallis in 1655, was 
also an infinite product: ^ 73 ^ 

7r /2 2 \ /4 4\ /6 6\ /8 8\ 

2 ~ vT ' 3 7 V3 ’ 5/ ‘ V5 ‘ 7/ V7 ’ 9/ 

The discovery of calculus, by English scientist Isaac Newton and German 
mathematician Gottfried Wilhelm Leibniz in the 1660s, led to the 
development of many infinite series for approximating 7T. Newton himself 
used an arcsin series to compute a 15 digit approximation of 7T in 1665 or 
1666, later writing "I am ashamed to tell you to how many figures I carried 
these computations, having no other business at the time."^ 72 ^ 



Isaac Newton used infinite 
series to compute n to 15 
digits, later writing "I am 
ashamed to tell you to how 
many figures I carried these 
computations".[ 72 1 


In Europe, Madhava's formula was rediscovered by Scottish mathematician James Gregory in 1671, and by 
Leibniz in 1674:^ 76 ^ 77 ^ 


arctan z = z — 


z 


7 



This formula, the Gregory-Leibniz series, equals n/4 when evaluated with Z = lP 7 ^ In 1699, English 
mathematician Abraham Sharp used the Gregory-Leibniz series for z = — to compute 7T to 71 digits, 

breaking the previous record of 39 digits, which was set with a polygonal algorithm. ^ 78 ^ The Gregory- 





































Leibniz for 2 = 1 series is simple, but converges very slowly (that is, approaches the answer gradually), so it 
is not used in modern 7T calculations.^ 


In 1706 John Machin used the Gregory-Leibniz series to produce an algorithm that converged much 
faster: [80] 


— = 4 arctan-arctan-. 

4 5 239 

Machin reached 100 digits of 71 with this formulaOther mathematicians created variants, now known as 
Machin-like formulae, that were used to set several successive records for calculating digits of 7lJ 81 ^ 
Machin-like formulae remained the best-known method for calculating 71 well into the age of computers, and 
were used to set records for 250 years, culminating in a 620-digit approximation in 1946 by Daniel 
Ferguson - the best approximation achieved without the aid of a calculating device.^ 


A record was set by the calculating prodigy Zacharias Dase, who in 1844 employed a Machin-like formula 
to calculate 200 decimals of 71 in his head at the behest of German mathematician Carl Friedrich GaussJ 83 ^ 
British mathematician William Shanks famously took 15 years to calculate 71 to 707 digits, but made a 
mistake in the 528th digit, rendering all subsequent digits incorrectJ 83 ^ 


Rate of convergence 

Some infinite series for 7T converge faster than others. Given the choice of two infinite series for 7T, 
mathematicians will generally use the one that converges more rapidly because faster convergence reduces 
the amount of computation needed to calculate 71 to any given accuracyJ 84 ^ A simple infinite series for 71 is 
the Gregory-Leibniz series:^ 



As individual terms of this infinite series are added to the sum, the total gradually gets closer to 71, and - 
with a sufficient number of terms - can get as close to n as desired. It converges quite slowly, though - after 
500,000 terms, it produces only five correct decimal digits of 7lJ 86 ^ 

An infinite series for 71 (published by Nilakantha in the 15th century) that converges more rapidly than the 

Gregory-Leibniz series is:{{sfn|Arndt|Haenel|2006|p=223|ps=: (formula 16.10). Note that (n - l)n(n + 1) = 
n 3_ n ;S7] 


7T = 3 + 


+ 


2x3x4 4x5x6 6x7x8 8x9x 10 

The following table compares the convergence rates of these two series: 


+ 


Infinite series for n 

After 

1st 

term 

After 

2nd 

term 

After 

3rd 

term 

After 

4th 

term 

After 

5th 

term 

Converges 

to: 

4 4 4 4 4 4 4 

1 3 + 5 7 + 9 U + 13 + " 

4.0000 

2.6666 

3.4666 

2.8952 

3.3396 

n - 3.1415 

4 4 4 

7F_ + 2x3x4 4x5x6 + 6x7x8 + 

3.0000 

3.1666 

3.1333 

3.1452 

3.1396 







After five terms, the sum of the Gregory-Leibniz series is within 0.2 of the correct value of 7T, whereas the 
sum of Nilakantha's series is within 0.002 of the correct value of 71. Nilakantha's series converges faster and 
is more useful for computing digits of 71. Series that converge even faster include Machin's series and 
Chudnovsky's series, the latter producing 14 correct decimal digits per term/ 84 ] 


Irrationality and transcendence 


Not all mathematical advances relating to 71 were aimed at increasing the accuracy of approximations. When 
Euler solved the Basel problem in 1735, finding the exact value of the sum of the reciprocal squares, he 
established a connection between 71 and the prime numbers that later contributed to the development and 
study of the Riemann zeta function/ 88 ] 



Swiss scientist Johann Heinrich Lambert in 1761 proved that 71 is irrational, meaning it is not equal to the 
quotient of any two whole numbers/ 19 ] Lambert's proof exploited a continued-fraction representation of the 
tangent function/ 89 ] French mathematician Adrien-Marie Legendre proved in 1794 that 7T 2 is also irrational. 
In 1882, German mathematician Ferdinand von Lindemann proved that 71 is transcendental, confirming a 
conjecture made by both Legendre and Euler/ 90 ^ 914 Hardy and Wright states that "the proofs were 
afterwards modified and simplified by Hilbert, Hurwitz, and other writers"/ 92 ] 


Adoption of the symbol n 


In the earliest usages, the Greek letter 71 was an abbreviation of the Greek 

word for periphery (nepupepeia)/ 93 ] and was combined in ratios with 5 (for 

diameter) or p (for radius) to form circle constantsJ 94 ^ 95 ^ 96 ^ (Before then, 

mathematicians sometimes used letters such as c or p instead/ 97 ]) The first 

recorded use is Oughtred's "5. 7r", to express the ratio of periphery and 

diameter in the 1647 and later editions of Clavis Mathematicae) 98 H 97 1 

Barrow likewise used "4" to represent the constant 3.14.../"] while 

0 

Gregory instead used to represent 6.28... .I 100 !! 95 ! 



Leonhard Euler popularized 
the use of the Greek letter n 
in works he published in 
1736 and 1748. 


The earliest known use of the Greek letter 71 alone to represent the ratio of a 
circle's circumference to its diameter was by Welsh mathematician William 
Jones in his 1706 work Synopsis Palmariorum Matheseos; or, a New 
Introduction to the Mathematics) 101 ^ { {sfn|Arndt|Haenel|2006|p=165|ps=: A 
facsimile of Jones' text is in Berggren, Borwein & Borwein 1997, pp. 108- 
109 The Greek letter first appears there in the phrase "1/2 Periphery (71)" in 
the discussion of a circle with radius one/ 102 ] However, he writes that his 
equations for 71 are from the "ready pen of the truly ingenious Mr. John 

Machin", leading to speculation that Machin may have employed the Greek letter before Jones/ 97 ] Jones' 
notation was not immediately adopted by other mathematicians, with the fraction notation still being used as 
late as 1767. [94][103 ] 


Euler started using the single-letter form beginning with his 1727 Essay Explaining the Properties of Air, 
though he used 71 = 6.28..., the ratio of radius to periphery, in this and some later writing/ 104 ^ 105 ] Euler first 
used 71 = 3.14... in his 1736 work Mechanica / 106 ] and continued in his widely-read 1748 work Introductio in 
analysin infinitorum (he wrote: "for the sake of brevity we will write this number as 71; thus 71 is equal to half 

























the circumference of a circle of radius I")/ 107 ] Because Euler corresponded heavily with other 
mathematicians in Europe, the use of the Greek letter spread rapidly, and the practice was universally 
adopted thereafter in the Western world/ 97 ] though the definition still varied between 3.14... and 6.28... as 
late as 176l/ 108 ] 

Modern quest for more digits 


Computer era and iterative algorithms 

The development of computers in the mid-20th century again revolutionized 
the hunt for digits of 7T. Mathematicians John Wrench and Levi Smith 
reached 1,120 digits in 1949 using a desk calculator/ 109 ] Using an inverse 
tangent (arctan) infinite series, a team led by George Reitwiesner and John 
von Neumann that same year achieved 2,037 digits with a calculation that 
took 70 hours of computer time on the ENIAC computer/ 110 ^ 111 ] The 
record, always relying on an arctan series, was broken repeatedly (7,480 
digits in 1957; 10,000 digits in 1958; 100,000 digits in 1961) until 1 million 
digits were reached in 1973/ 110 ] 

Two additional developments around 1980 once again accelerated the ability 
to compute 71. First, the discovery of new iterative algorithms for computing 
71, which were much faster than the infinite series; and second, the invention 
of fast multiplication algorithms that could multiply large numbers very 
rapidly/ 112 ] Such algorithms are particularly important in modern 7T 
computations because most of the computer's time is devoted to 
multiplication/ 113 ] They include the Karatsuba algorithm, 

Toom-Cook multiplication, and Fourier transform-based 
methods / 114 ] 

The iterative algorithms were independently published in 1975- 
1976 by physicist Eugene Salamin and scientist Richard 
Brent/ 115 ] These avoid reliance on infinite series. An iterative 
algorithm repeats a specific calculation, each iteration using the 
outputs from prior steps as its inputs, and produces a result in 
each step that converges to the desired value. The approach was 
actually invented over 160 years earlier by Carl Friedrich 
Gauss, in what is now termed the arithmetic-geometric mean 
method (AGM method) or Gauss-Legendre algorithm/ 115 ] As 
modified by Salamin and Brent, it is also referred to as the 
Brent-Salamin algorithm. 

The iterative algorithms were widely used after 1980 because 
they are faster than infinite series algorithms: whereas infinite series typically increase the number of correct 
digits additively in successive terms, iterative algorithms generally multiply the number of correct digits at 
each step. For example, the Brent-Salamin algorithm doubles the number of digits in each iteration. In 1984, 
brothers John and Peter Borwein produced an iterative algorithm that quadruples the number of digits in 
each step; and in 1987, one that increases the number of digits five times in each step/ 116 ] Iterative methods 
were used by Japanese mathematician Yasumasa Kanada to set several records for computing 7T between 
1995 and 2002/ liy ] This rapid convergence comes at a price: the iterative algorithms require significantly 
more memory than infinite series/ 117 ] 


The Gauss-Legendre iterative algorithm: 
Initialize 

a o=l, 6 ° = yf’ to= J’ Po =1 - 

Iterate 

_ a n+^n u _ /t—— 

a n+ 1— 2 ’ ^n+1— 

^n+l = ^n~ Pni a n~ a n+l) ? Pn+l =( ^Pn‘ 

Then an estimate for 7T is given by 

(an+bn)^ 1 

7TRi---. 

4 tn 



John von Neumann was 
part of the team that first 
used a digital computer, 
ENIAC, to compute n. 




Motives for computing n 


For most numerical calculations 
involving 71, a handful of digits provide 
sufficient precision. According to Jorg 
Arndt and Christoph Haenel, thirty-nine 
digits are sufficient to perform most 
cosmological calculations, because that 
is the accuracy necessary to calculate the 
circumference of the observable universe 
with 


Record approximations of pf 
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As mathematicians discovered new algorithms, and computers 
became available, the number of known decimal digits of n 
increased dramatically. The vertical scale is logarithmic. 


a precision of one atom. L 
Accounting for additional digits needed 
to compensate for computational round¬ 
off errors, Arndt concludes that a few 
hundred digits would suffice for any 
scientific application. Despite this, 
people have worked strenuously to 
compute 7T to thousands and millions of 

digits/ 119 ^ This effort may be partly ascribed to the human compulsion to break records, and such 
achievements with 7T often make headlines around the world/ 120 -^ 121 ! They also have practical benefits, such 
as testing supercomputers, testing numerical analysis algorithms (including high-precision multiplication 
algorithms); and within pure mathematics itself, providing data for evaluating the randomness of the digits 

Of 7T.t 122 ] 


Rapidly convergent series 


Modern 71 calculators do not use iterative algorithms exclusively. New 
infinite series were discovered in the 1980s and 1990s that are as fast as 
iterative algorithms, yet are simpler and less memory intensive/ 117 ! The fast 
iterative algorithms were anticipated in 1914, when the Indian 
mathematician Srinivasa Ramanujan published dozens of innovative new 
formulae for 71, remarkable for their elegance, mathematical depth, and rapid 
convergence/ 123 ! One of his formulae, based on modular equations, is 


1 

7T 


2^2 

9801 


E 

k =0 


(4fc)!(1103 + 26390&) 
A:! 4 (396 4fc ) 


This series converges much more rapidly than most arctan series, including 
Machin's formula/ 124 ! Bill Gosper was the first to use it for advances in the 
calculation of 7T, setting a record of 17 million digits in 1985/ 125 ! 
Ramanujan's formulae anticipated the modern algorithms developed by the 
Borwein brothers and the Chudnovsky brothers/ 126 ! The Chudnovsky 
formula developed in 1987 is 



Srinivasa Ramanujan, 
working in isolation in India, 
produced many innovative 
series for computing n. 


1 

7T 


12 

640320 3 / 2 


OO 


E 


(6fc) 1(13591409 + 545140134&) 
(3/c)!(/c!) 3 (—640320) 3fc 


It produces about 14 digits of n per term/ 127 ! anc j h as b een usec j f or several record-setting 7T calculations, 
including the first to surpass 1 billion (10 9 ) digits in 1989 by the Chudnovsky brothers, 2.7 trillion 
(2.7xl0 12 ) digits by Fabrice Bellard in 2009/ 128 ! 10 trillion (10 13 ) digits in 2011 by Alexander Yee and 





















Shigeru Rondo/ 129 ! anc [ over 22 trillion digits in 2016 by Peter Trueb/ 130 ^ 131 ! For similar formulas, see 
also the Ramanujan-Sato series. 


In 2006, mathematician Simon Plouffe used the PSLQ integer relation algorithm/ 132 ! to generate several new 
formulas for 7T, conforming to the following template: 

^ f 1/ ° , b _ , c 

“j n k \Q n ~ 1 q 2n — 1 g 4n — 1 

where q is e n (Gelfond's constant), k is an odd number, and a, b, C are certain rational numbers that Plouffe 
computed/ 133 ! 



Monte Carlo methods 


Monte Carlo methods, which 
evaluate the results of multiple 
random trials, can be used to 
create approximations of 
7i/ 134 ! Buffon's needle is one 
such technique: If a needle of 
length £ is dropped n times on 
a surface on which parallel 
lines are drawn f units apart, 
and if X of those times it comes 
to rest crossing a line (X > 0), 
then one may approximate 7T 
based on the counts/ 135 ! 





Buffon's needle. Needles a 
and b are dropped 
randomly. 


Random dots 
are placed on 
the quadrant of 
a square with a 
circle inscribed 
in it. 


Monte Carlo methods, based on random trials, can be used to approximate n. 


2 ni 

7r ~ -. 

xt 

Another Monte Carlo method for computing 7T is to draw a circle inscribed in a square, and randomly place 
dots in the square. The ratio of dots inside the circle to the total number of dots will approximately equal 

ji/4 [136] 


Another way to calculate 71 using probability is to start with a 
random walk, generated by a sequence of (fair) coin tosses: 
independent random variables X^ such that £ {-1,1} with equal 
probabilities. The associated random walk is 

n 

w n = Y,x t 

k =1 

so that, for each n, W n is drawn from a shifted and scaled binomial 
distribution. As n varies, W n defines a (discrete) stochastic process. 
Then 7T can be calculated by^ 137 ! 



Five random walks with 200 steps. 
The sample mean of | W 200 I' s 
p = 56/5, and so 2(200)p“ 2 « 3.19 
is within 0.05 of n. 


7T = lim 

n—00 



















This Monte Carlo method is independent of any relation to circles, and is a consequence of the central limit 
theorem, discussed below. 

These Monte Carlo methods for approximating 7T are very slow compared to other methods, and do not 
provide any information on the exact number of digits that are obtained. Thus they are never used to 
approximate 71 when speed or accuracy is desired/ 138 ] 


Spigot algorithms 

Two algorithms were discovered in 1995 that opened up new avenues of research into 71. They are called 
spigot algorithms because, like water dripping from a spigot , they produce single digits of 71 that are not 
reused after they are calculated/ 139 ** 440 ] This is in contrast to infinite series or iterative algorithms, which 
retain and use all intermediate digits until the final result is produced/ 139 * 

Mathematicians Stan Wagon and Stanley Rabinowitz produced a simple spigot algorithm in 
1995/ 140 ^ 141 ^ 142 ] Its speed is comparable to arctan algorithms, but not as fast as iterative algorithms/ 141 * 

Another spigot algorithm, the BBP digit extraction algorithm, was discovered in 1995 by Simon 
Plouffe/ 143 **^ 44 * 


_ _ V 1 ( 4 2 1 1 \ 

^ _ ^ 16* \ 8A; + 1 ~~ 8fc + 4 ~~ 8fc + 5 ~~ 8k + 6 ) ' 

This formula, unlike others before it, can produce any individual hexadecimal digit of 71 without calculating 
all the preceding digits/ 143 -* Individual binary digits may be extracted from individual hexadecimal digits, 
and octal digits can be extracted from one or two hexadecimal digits. Variations of the algorithm have been 
discovered, but no digit extraction algorithm has yet been found that rapidly produces decimal digits/ 145 -* An 
important application of digit extraction algorithms is to validate new claims of record 71 computations: After 
a new record is claimed, the decimal result is converted to hexadecimal, and then a digit extraction algorithm 
is used to calculate several random hexadecimal digits near the end; if they match, this provides a measure of 
confidence that the entire computation is correct/ 129 * 

Between 1998 and 2000, the distributed computing project PiHex used Bellard's formula (a modification of 
the BBP algorithm) to compute the quadrillionth (10 15 th) bit of 71, which turned out to be O.* 146 * In 
September 2010, a Yahoo! employee used the company's Hadoop application on one thousand computers 
over a 23-day period to compute 256 bits of 71 at the two-quadrillionth (2xl0 15 th) bit, which also happens to 
be zero/ 147 ] 

Role and characterizations in mathematics 


Because 7T is closely related to the circle, it is found in many formulae from the fields of geometry and 
trigonometry, particularly those concerning circles, spheres, or ellipses. Other branches of science, such as 
statistics, physics, Fourier analysis, and number theory, also include 71 in some of their important formulae. 


Geometry and trigonometry 

7i appears in formulae for areas and volumes of geometrical shapes based on circles, such as ellipses, 
spheres, cones, and tori. Below are some of the more common formulae that involve Tl/ 148 * 

■ The circumference of a circle with radius r is 2 nr. 





■ The area of a circle with radius r is nr 2 . 

4 ^ 

■ The volume of a sphere with radius r is ^nr . 

■ The surface area of a sphere with radius r is 4nr 2 . 

The formulae above are special cases of the volume of the n- 
dimensional ball and the surface area of its boundary, the (n-1)- 
dimensional sphere, given below. 

Definite integrals that describe circumference, area, or volume of 
shapes generated by circles typically have values that involve 7T. For 
example, an integral that specifies half the area of a circle of radius 
one is given by:^ 149 ^ 

f ^/l — x 2 dx = —. 

J -1 2 



The area of the circle equals n times 
the shaded area. 


In that integral the function Vl - X 2 represents the top half of a circle (the square root is a consequence of 

1 

the Pythagorean theorem), and the integral f computes the area between that half of a circle and the x 
axis. 


The trigonometric functions rely on angles, and 
mathematicians generally use radians as units of 
measurement. 7T plays an important role in angles 
measured in radians, which are defined so that a 
complete circle spans an angle of 271 radians.^ 150 ^ 
The angle measure of 180° is equal to n radians, 
and 1° = 71/180 radiansJ 150 ^ 


Common trigonometric functions have periods 
that are multiples of 71; for example, sine and 
cosine have period 27r/ 151 ^ so for any angle 9 and 
any integer k, 



Sine and cosine functions repeat with period 2n. 


sin# = sin(# + 2irk) and cos 9 = cos (6 + 27rfc)J 151 ^ 


Eigenvalues 

Many of the appearances of n in the formulas of mathematics and the sciences have to do with its close 
relationship with geometry. However, 71 also appears in many natural situations having apparently nothing to 
do with geometry. 

In many applications, it plays a distinguished role as an eigenvalue. For example, an idealized vibrating 
string can be modelled as the graph of a function f on the unit interval [0,1], with fixed ends 
f(0) = f(l) = 0. The modes of vibration of the string are solutions of the differential equation 
f"{x) + \f{x) — 0, or f" (t) — —A f(x). Thus A is an eigenvalue of the second derivative operator 
/ i —y f", and is constrained by Sturm-Liouville theory to take on only certain specific values. It must be 

positive, since the operator is negative definite, so it is convenient to write A = v , where v > 0 is called the 
wavenumber. Then f(x ) = sin(Ti x) satisfies the boundary conditions and the differential equation with 
v — nJ 152 ] 



























The value 71 is, in fact, the least such value of the wavenumber, and 
is associated with the fundamental mode of vibration of the string. 
One way to show this is by estimating the energy, which satisfies 
Wirtinger's inequality^ 153 ] for a function f : [0, 1] -> C with 
f(0) = f( 1) = 0 and f, f' both square integrable, we have: 



f 


i r / \f(x)\ dx < 




with equality precisely when f is a multiple of sin(TT v). Here n 
appears as an optimal constant in Wirtinger's inequality, and it 
follows that it is the smallest wavenumber, using the variational 
characterization of the eigenvalue. As a consequence, 7T is the 
smallest singular value of the derivative operator on the space of 
functions on [0,1] vanishing at both endpoints (the Sobolev space 

H] [ 0 . 1 ]). 



The overtones of a vibrating string 
are eigenfunctions of the second 
derivative, and form a harmonic 
progression. The associated 
eigenvalues form the arithmetic 
progression of integer multiples of n. 


Inequalities 


The number 71 serves appears in similar eigenvalue problems in 
higher-dimensional analysis. As mentioned above, it can be 
characterized via its role as the best constant in the isoperimetric 
inequality : the area A enclosed by a plane Jordan curve of perimeter 
P satisfies the inequality 

4t tA < P 2 , 

and equality is clearly achieved for the circle, since in that case 
A = nr 2 and P = 2nr. [154] 

Ultimately as a consequence of the isoperimetric inequality, 7T 
appears in the optimal constant for the critical Sobolev inequality in 
n dimensions, which thus characterizes the role of 71 in many 
physical phenomena as well, for example those of classical potential 
theoryj 155 ^ 156 ^ 157 ] In two dimensions, the critical Sobolev 
inequality is 

2*11/112 < IIV/IU 

for f a smooth function with compact support in R , V/ is the 
gradient of f, and ||/||2 and |]V/||i refer respectively to the L and 



The ancient city of Carthage was 
the solution to an isoperimetric 
problem, according to a legend 
recounted by Lord Kelvin 
(Thompson 1894): those lands 
bordering the sea that Queen Dido 
could enclose on all other sides 
within a single given oxhide, cut into 
strips. 


4 

L -norm. The Sobolev inequality is equivalent to the isoperimetric inequality (in any dimension), with the 
same best constants. 


Wirtinger's inequality also generalizes to higher-dimensional Poincare inequalities that provide best 
constants for the Dirichlet energy of an n-dimensional membrane. Specifically, 7T is the greatest constant 
such that 



for all convex subsets G of R n of diameter 1, and square-integrable functions u on G of mean zero, f 158 ^ Just 
as Wirtinger's inequality is the variational form of the Dirichlet eigenvalue problem in one dimension, the 
Poincare inequality is the variational form of the Neumann eigenvalue problem, in any dimension. 


Fourier transform and Heisenberg uncertainty principle 


The constant 7T also appears as a critical spectral parameter in the 
Fourier transform. This is the integral transform, that takes a 
complex-valued integrable function f on the real line to the function 
defined as: 

/ OO 

f(x)e~ 2 ™( dx. 

-oo 

Although there are several different conventions for the Fourier 
transform and its inverse, any such convention must involve 7T 
somewhere. The above is the most canonical definition, however, 

giving the unique unitary operator on L that is also an algebra 
homomorphism of L * to L 00 .^ 59 ! 

The Heisenberg uncertainty principle also contains the number 7T. 
The uncertainty principle gives a sharp lower bound on the extent to 
which it is possible to localize a function both in space and in 
frequency: with our conventions for the Fourier transform. 



An animation of a geodesic in the 
Heisenberg group, showing the 
close connection between the 
Heisenberg group, isoperimetry, and 
the constant n. The cumulative 
height of the geodesic is equal to 
the area of the shaded portion of the 
unit circle, while the arc length (in 
the Carnot-Caratheodory metric) is 
equal to the circumference. 


(J^x 2 \ f ( X )\ 2 dx') (£Vi/(oi 2 d?) > (^/jAz)i 2 <&) 2 . 

The physical consequence, about the uncertainty in simultaneous position and momentum observations of a 
quantum mechanical system, is discussed below. The appearance of 71 in the formulae of Fourier analysis is 
ultimately a consequence of the Stone-von Neumann theorem, asserting the uniqueness of the Schrodinger 
representation of the Heisenberg group J 160 ^ 


Gaussian integrals 



















The fields of probability and statistics frequently use the normal 
distribution as a simple model for complex phenomena; for example, 
scientists generally assume that the observational error in most 
experiments follows a normal distribution.^ 161 ^ The Gaussian 
function, which is the probability density function of the normal 
distribution with mean |i and standard deviation cf, naturally contains 

Jf;[ 162 ] 


/(*) 


1 e -( X ~V) 2 /( 2 °' 2 ) ' 

ay/ 27T 


The factor of —p= makes the area under the graph of f equal to one, 

Y 27T 

as is required for a probability distribution. This follows from a 
change of variables in the Gaussian integral^ 1621 



-2 -1 □ 1 2 


A graph of the Gaussian function 

_ 2 

f(x) = e x . The coloured region 
between th£function and the x-axis 
has area Vji. 



which says that the area under the basic bell curve in the figure is equal to the square root of 7T. 


The central limit theorem explains the central role of normal 
distributions, and thus of 7T, in probability and statistics. This 
theorem is ultimately connected with the spectral characterization of 
71 as the eigenvalue associated with the Heisenberg uncertainty 
principle, and the fact that equality holds in the uncertainty principle 
only for the Gaussian functionJ 163 ^ Equivalently, 7T is the unique 

TTX 2 

constant making the Gaussian normal distribution e equal to its 
own Fourier transform.^ 164 ! Indeed, according to Howe (1980), the 
"whole business" of establishing the fundamental theorems of 
Fourier analysis reduces to the Gaussian integral. 



7i can be computed from the 
distribution of zeros of a one¬ 
dimensional Wiener process 


Projective geometry 

Let V be the set of all twice differentiable real functions / : R —)■ R that satisfy the ordinary differential 
equation f" (x) + f(x) = 0. Then V is a two-dimensional real vector space, with two parameters 
corresponding to a pair of initial conditions for the differential equation. For any teR, let e* : V —> R be 
the evaluation functional, which associates to each / e V the value e*(/) = f(t) of the function f at the 
real point t. Then, for each t, the kernel of e* is a one-dimensional linear subspace of V. Hence t H* ker e t 
defines a function from M —> P(V) from the real line to the real projective line. This function is periodic, 
and the quantity 71 can be characterized as the period of this mapJ 165 ^ 


Topology 


The constant 7T appears in the Gauss-Bonnet formula which relates the differential geometry of surfaces to 
their topology . Specifically, if a compact surface Z has Gauss curvature K, then 



KdA = 2ttx(S) 























where/(X) is the Euler characteristic, which is an integer.^ 166 ! An 
example is the surface area of a sphere S of curvature 1 (so that its 
radius of curvature, which coincides with its radius, is also 1.) The 
Euler characteristic of a sphere can be computed from its homology 
groups and is found to be equal to two. Thus we have 

A(S ) — f 1 dA — • 2 = 47t 

Js 

reproducing the formula for the surface area of a sphere of radius 1. 

The constant appears in many other integral formulae in topology, in 
particular, those involving characteristic classes via the Chern-Weil 
homomorphism^ 167 ^ 

Vector calculus 

Vector calculus is a branch of calculus that is concerned with the 
properties of vector fields, and has many physical applications such 
as to electricity and magnetism. The Newtonian potential for a point 
source Q situated at the origin of a three-dimensional Cartesian 
coordinate system is'- 168 ^ 



Uniformization of the Klein quartic, a 
surface of genus three and Euler 
characteristic -4, as a quotient of 
the hyperbolic plane by the 
symmetry group PSL(2,7) of the 
Fano plane. The hyperbolic area of 
a fundamental domain is 8 ti, by 
Gauss-Bonnet. 


V(x) 


kQ 


x 


which represents the potential energy of a unit mass (or charge) 
placed a distance |x| from the source, and k is a dimensional 
constant. The field, denoted here by E, which may be the 
(Newtonian) gravitational field or the (Coulomb) electric field, is the 
negative gradient of the potential: 

E = - W. 



The techniques of vector calculus 
can be understood in terms of 
decompositions into spherical 
harmonics (shown) 


Special cases include Coulomb's law and Newton's law of universal 
gravitation. Gauss' law states that the outward flux of the field 

through any smooth, simple, closed, orientable surface S containing the origin is equal to 4nkQ: 


4irkQ = (If) E • dA. 
JJs 


It is standard to absorb this factor of 4 tt into the constant k, but this argument shows why it must appear 
somewhere. Furthermore, 4tt is the surface area of the unit sphere, but we have not assumed that S is the 
sphere. However, as a consequence of the divergence theorem, because the region away from the origin is 

vacuum (source-free) it is only the homology class of the surface S in R J \{0} that matters in computing the 
integral, so it can be replaced by any convenient surface in the same homology class, in particular, a sphere, 
where spherical coordinates can be used to calculate the integral. 


A consequence of the Gauss law is that the negative Laplacian of the potential V is equal to 4nkQ times the 
Dirac delta function: 


AE(x) = —4:7rkQ8(x.). 





























More general distributions of matter (or charge) are obtained from this by convolution, giving the Poisson 
equation 

AV(x) = —47 t fcp(x) 
where p is the distribution function. 

The constant 7T also plays an analogous role in four-dimensional 
potentials associated with Einstein's equations, a fundamental 
formula which forms the basis of the general theory of relativity and 
describes the fundamental interaction of gravitation as a result of 
spacetime being curved by matter and energy : ^ 169 ^ 

1 87 tG 

■ti-fw 1Z ^9/u' AfiV — 7 ^ t u/ ’ 

A C 4 

where R„ v is the Ricci curvature tensor, R is the scalar curvature, 
g^y is the metric tensor, A is the cosmological constant, G is Newton's gravitational constant, C is the speed 
of light in vacuum, and is the stress-energy tensor. The left-hand side of Einstein's equation is a non¬ 
linear analogue of the Laplacian of the metric tensor, and reduces to that in the weak field limit, with the A g 
term playing the role of a Lagrange multiplier , and the right-hand side is the analogue of the distribution 
function, times 8 tt. 



Einstein's equation states that the 
curvature of space-time is produced 
by the matter-energy content. 


Cauchy's integral formula 


One of the key tools in complex analysis is contour integration of a 
function over a positively oriented (rectifiable) Jordan curve y. A 
form of Cauchy's integral formula states that if a point Zn is interior 

to Y , then [170] 


z- z 0 

Although the curve y is not a circle, and hence does not have any 
obvious connection to the constant 7T, a standard proof of this result 
uses Morera's theorem, which implies that the integral is invariant 
under homotopy of the curve, so that it can be deformed to a circle 
and then integrated explicitly in polar coordinates. More generally, it 
is true that if a rectifiable closed curve y does not contain Zq, then 
the above integral is 2ni times the winding number of the curve. 
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Complex analytic functions can be 
visualized as a collection of 
streamlines and equipotentials, 
systems of curves intersecting at 
right angles. Here illustrated is the 
complex logarithm of the Gamma 
function. 


The general form of Cauchy's integral formula establishes the 

relationship between the values of a complex analytic function f(z) on the Jordan curve y and the value of 
f(z) at any interior point Zq of y;[ 171 ][ 172 l 


- dz = 27rif(zo) 

Z~ Zq 

provided f(z) is analytic in the region enclosed by y and extends continuously to y. Cauchy's integral 
formula is a special case of the residue theorem, that if g(z ) is a meromorphic function the region enclosed 
by y and is continuous in a neighbourhood of y, then 































(P g(z) dz = 2ni^2Res(g, a/-) 

J 7 


where the sum is of the residues at the poles of g(z). 


The gamma function and Stirling's approximation 


The factorial function n ! is the product of all of the positive integers 
through n. The gamma function extends the concept of factorial 
(normally defined only for non-negative integers) to all complex 
numbers, except the negative real integers. When the gamma 
function is evaluated at half-integers, the result contains 71; for 

example T(l/2) = yTr and r(5/2) = ^^.1 173 1 


The gamma function is defined by its Weierstrass product 
development! 174 ] 


r(*) 



e z / n 

1 + zjn 


where y is the Euler-Mascheroni constant. Evaluated at Z - 1/2 and 

squared, the equation T(l/2) z = tt reduces to the Wallis product 
formula. The gamma function is also connected to the Riemann zeta 
function and identities for the functional determinant, in which the 
constant 7T plays an important role. 



The Hopf fibration of the 3-sphere, 
by Villarceau circles, over the 
complex projective line with its 
Fubini-Study metric (three parallels 
are shown). The identity 
S 3 (l)/S 2 (l) = ti/ 2 is a 
consequence. 


The gamma function is used to calculate the volume V n (r ) of the n-dimensional ball of radius r in Euclidean 
n-dimensional space, and the surface area S n -l( r ) °f its boundary, the (n-l)-dimensional sphere! 175 ] 


Vn(r) = 


7 r 


.71 / 2 


>n -1 


m = 


r(! + i) ’ 

.71 / 2 


n tv 


n— 1 


r(f+ i) 


Further, it follows from the functional equation that 


27rr = 


S„+i(r) 
KM ' 


The gamma function can be used to create a simple approximation to the factorial function n! for large n : 
ri! ~ \/2ivn (^) ” which is known as Stirling's approximation! 176 ] Equivalently, 


7r = lim 


e n\ 


ti^oo 2iTl^ n ^ 


















As a geometrical application of Stirling's approximation, let A n denote the standard simplex in n- 
dimensional Euclidean space, and (n + l)A n denote the simplex having all of its sides scaled up by a factor 
of n + 1. Then 


Vol((n + l)A n ) = 


(n + 1 ) r 


n! 


0 n +1 


V27m 


Ehrhart's volume conjecture is that this is the (optimal) upper bound on the volume of a convex body 
containing only one lattice point J 177 -^ 


Number theory and Riemann zeta function 


The Riemann zeta function £(s) is used in many areas of 
mathematics. When evaluated at S = 2 it can be written as 



Finding a simple solution for this infinite series was a famous 
problem in mathematics called the Basel problem. Leonhard Euler 

solved it in 1735 when he showed it was equal to TI 2 /6J 88 ^ Euler's 
result leads to the number theory result that the probability of two 
random numbers being relatively prime (that is, having no shared 

factors) is equal to 6 /ti 2 J 178 ^ 179 ^ This probability is based on the 
observation that the probability that any number is divisible by a 
prime p is 1/p (for example, every 7th integer is divisible by 7.) 
Hence the probability that two numbers are both divisible by this 

prime is 1/p , and the probability that at least one of them is not is 

1 ~ 1/p . For distinct primes, these divisibility events are mutually 
independent; so the probability that two numbers are relatively prime 
is given by a product over all primes :^ 180 ^ 
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Each prime has an associated 
Priifer group, which are arithmetic 
localizations of the circle. The L- 
functions of analytic number theory 
are also localized in each prime p. 



Solution of the Basel problem using 
the Weil conjecture: the value of 
<"(2) is the hyperbolic area of a 
fundamental domain of the modular 
group, times 7.n. 


This probability can be used in conjunction with a random number generator to approximate 71 using a 
Monte Carlo approachJ 181 ^ 

The solution to the Basel problem implies that the geometrically derived quantity 7T is connected in a deep 
way to the distribution of prime numbers. This is a special case of Weil's conjecture on Tamagawa numbers, 
which asserts the equality of similar such infinite products of arithmetic quantities, localized at each prime p, 
and a geometrical quantity: the reciprocal of the volume of a certain locally symmetric space. In the case of 
the Basel problem, it is the hyperbolic 3-manifold SL, 2 (R)/SL, 2 (Z)J 182 ^ 





























The zeta function also satisfies Riemann's functional equation, which involves 71 as well as the gamma 
function: 

C(s) = 2 s 7T s-1 s in (y) T(1 - s) C(1 - s). 

Furthermore, the derivative of the zeta function satisfies 
exp(-C , (0)) = y/2ir. 

A consequence is that 71 can be obtained from the functional determinant of the harmonic oscillator. This 
functional determinant can be computed via a product expansion, and is equivalent to the Wallis product 
formula.^ 183 ! The calculation can be recast in quantum mechanics, specifically the variational approach to 
the spectrum of the hydrogen atomJ 184 ^ 


Fourier series 


The constant 7T also appears naturally in Fourier series of periodic 
functions. Periodic functions are functions on the group T =R/Z of 
fractional parts of real numbers. The Fourier decomposition shows 
that a complex-valued function f on T can be written as an infinite 
linear superposition of unitary characters of T. That is, continuous 
group homomorphisms from T to the circle group 1/(1) of unit 
modulus complex numbers. It is a theorem that every character of T 
is one of the complex exponentials e n (x) = e 2mnx . 

There is a unique character on T, up to complex conjugation, that is 
a group isomorphism. Using the Haar measure on the circle group, 
the constant 7T is half the magnitude of the Radon-Nikodym 
derivative of this character. The other characters have derivatives 
whose magnitudes are positive integral multiples of 271J 18 ^ As a 
result, the constant n is the unique number such that the group T, 
equipped with its Haar measure, is Pontrjagin dual to the lattice of 
integral multiples of 271J 186 ^ This is a version of the one-dimensional 
Poisson summation formula. 



n appears in characters of p-adic 
numbers (shown), which are 
elements of a Priifer group. Tate's 
thesis makes heavy use of this 
machinery.! 185 ^ 


Modular forms and theta functions 


The constant 71 is connected in a deep way with the theory of 
modular forms and theta functions. For example, the Chudnovsky 
algorithm involves in an essential way the j-invariant of an elliptic 
curve. 

Modular forms are holomorphic functions in the upper half plane 
characterized by their transformation properties under the modular 
group SL 2 (Z) (or its various subgroups), a lattice in the group 
SL 2 (M). An example is the Jacobi theta function 

00 

^ ^ g2ninz+iirn? t 
n=— 00 



Theta functions transform under the 
lattice of periods of an elliptic curve. 


























which is a kind of modular form called a Jacobi formJ 187 ^ This is sometimes written in terms of the nome 


The constant 71 is the unique constant making the Jacobi theta function an automorphic form, which means 
that it transforms in a specific way. Certain identities hold for all automorphic forms. An example is 

0{z + t,t) = e- niT - 2niz 6(zi t), 

which implies that 0 transforms as a representation under the discrete Heisenberg group . General modular 
forms and other theta functions also involve 71, once again because of the Stone-von Neumann theorem.t 187 ] 


Cauchy distribution and potential theory 


The Cauchy distribution 


g( x ) 


i i 

7T X 2 + 1 


is a probability density function. The total probability is equal to one, 
owing to the integral: 




The Witch of Agnesi, named for 
Maria Agnesi (1718-1799), is a 
geometrical construction of the 
graph of the Cauchy distribution. 


The Shannon entropy of the Cauchy distribution is equal to ln(4TT), 
which also involves 7T. 


The Cauchy distribution plays an important role in potential theory because it is the 
simplest Furstenberg measure, the classical Poisson kernel associated with a 
Brownian motion in a half-planeJ 188 ^ Conjugate harmonic functions and so also 
the Hilbert transform are associated with the asymptotics of the Poisson kernel. The 
Hilbert transform H is the integral transform given by the Cauchy principal value 
of the singular integral 



fjx ) dx 
x — t 


The constant 7T is the unique (positive) normalizing factor such that H defines a 
linear complex structure on the Hilbert space of square-integrable real-valued 
functions on the real lineJ 189 ^ The Hilbert transform, like the Fourier transform, 
can be characterized purely in terms of its transformation properties on the Hilbert 

space L (R) : up to a normalization factor, it is the unique bounded linear operator 
that commutes with positive dilations and anti-commutes with all reflections of 
constant 7T is the unique normalizing factor that makes this transformation unitary. 



The Cauchy 
distribution governs 
the passage of 
Brownian particles 
through a 
membrane. 


the real lineJ 190 ^ The 


Complex dynamics 

An occurrence of 7T in the Mandelbrot set fractal was discovered by David Boll in 199lJ 191 ^ He examined 
the behaviour of the Mandelbrot set near the "neck" at (-0.75, 0). If points with coordinates (-0.75, 8) are 
considered, as 8 tends to zero, the number of iterations until divergence for the point multiplied by 8 




























converges to 71. The point ( 0.25 + £, 0 ) at the cusp of the large 
"valley" on the right side of the Mandelbrot set behaves similarly: 
the number of iterations until divergence multiplied by the square 
root of £ tends to 7 r.[ 191 l[ 192 ] 

Outside mathematics 


Describing physical phenomena 

Although not a physical constant, 71 appears routinely in equations 
describing fundamental principles of the universe, often because of 
7T's relationship to the circle and to spherical coordinate systems. A 
simple formula from the field of classical mechanics gives the 
approximate period T of a simple pendulum of length L, swinging 
with a small amplitude (g is the earth's gravitational acceleration):^ 19 ^ 

T « 2tt 

One of the key formulae of quantum mechanics is Heisenberg's uncertainty principle, which shows that the 
uncertainty in the measurement of a particle's position (Ax) and momentum (Ap) cannot both be arbitrarily 
small at the same time (where h is Planck's constant)^ 194 ! 

. . h 

Ax A p > —. 

47r 

The fact that 71 is approximately equal to 3 plays a role in the relatively long lifetime of orthopositronium. 
The inverse lifetime to lowest order in the fine-structure constant a is^ 195 ^ 



where m is the mass of the electron. 

71 is present in some structural engineering formulae, such as the buckling formula derived by Euler, which 
gives the maximum axial load F that a long, slender column of length L, modulus of elasticity E, and area 
moment of inertia I can carry without buckling: ^ 196 ^ 



The field of fluid dynamics contains 7T in Stokes' law, which approximates the frictional force F exerted on 
small, spherical objects of radius R, moving with velocity V in a fluid with dynamic viscosity 

F = Girr/Rv. 

In electromagnetics, the vacuum permeability constant p 0 appears in Maxwell's equations, which describe 
the properties of electric and magnetic fields and electromagnetic radiation. Before 20 May 2019, it was 
defined as exactly 




7i can be computed from the 
Mandelbrot set, by counting the 
number of iterations required before 
point (-0.75, e) diverges. 


























Ho = 4tt x 10“ 7 H/m « 1.2566370614... x 10“ 6 N/A 2 . 


A relation for the speed of light in vacuum, C can be derived from Maxwell's equations in the medium of 
classical vacuum using a relationship between p 0 an d the electric constant (vacuum permittivity), £q in SI 
units: 


1 

c = . 

Under ideal conditions (uniform gentle slope on a homogeneously erodible substrate), the sinuosity of a 
meandering river approaches 71. The sinuosity is the ratio between the actual length and the straight-line 
distance from source to mouth. Faster currents along the outside edges of a river's bends cause more erosion 
than along the inside edges, thus pushing the bends even farther out, and increasing the overall loopiness of 
the river. However, that loopiness eventually causes the river to double back on itself in places and "short- 
circuit", creating an ox-bow lake in the process. The balance between these two opposing factors leads to an 
average ratio of 71 between the actual length and the direct distance between source and mouth/ 198 ^ 1 "] 


Memorizing digits 

Piphilology is the practice of memorizing large numbers of digits of 7t/ 20 °] and world-records are kept by 
the Guinness World Records. The record for memorizing digits of 71, certified by Guinness World Records, is 
70,000 digits, recited in India by Rajveer Meena in 9 hours and 27 minutes on 21 March 2015/ 201 ] In 2006, 
Akira Haraguchi, a retired Japanese engineer, claimed to have recited 100,000 decimal places, but the claim 
was not verified by Guinness World Records/ 202 ] 

One common technique is to memorize a story or poem in which the word lengths represent the digits of 7T: 
The first word has three letters, the second word has one, the third has four, the fourth has one, the fifth has 
five, and so on. Such memorization aids are called mnemonics. An early example of a mnemonic for pi, 
originally devised by English scientist James Jeans, is "How I want a drink, alcoholic of course, after the 
heavy lectures involving quantum mechanics."1 20 °1 When a poem is used, it is sometimes referred to as a 
piem. Poems for memorizing 7T have been composed in several languages in addition to English/ 200 ] 
Record-setting 71 memorizers typically do not rely on poems, but instead use methods such as remembering 
number patterns and the method of loci/ 203 ] 

A few authors have used the digits of n to establish a new form of constrained writing, where the word 
lengths are required to represent the digits of 71. The Cadaeic Cadenza contains the first 3835 digits of 71 in 

this manner/ 204 ^ and the full-length book Not a Wake contains 10,000 words, each representing one digit of 
^ [205] 

In popular culture 

Perhaps because of the simplicity of its definition and its ubiquitous presence in formulae, 71 has been 
represented in popular culture more than other mathematical constructs/ 206 ] 

In the 2008 Open University and BBC documentary co-production, The Story of Maths, aired in October 
2008 on BBC Four, British mathematician Marcus du Sautoy shows a visualization of the - historically first 
exact - formula for calculating 71 when visiting India and exploring its contributions to trigonometry/ 207 ] 

In the Palais de la Decouverte (a science museum in Paris) there is a circular room known as the pi room. On 
its wall are inscribed 707 digits of 71. The digits are large wooden characters attached to the dome-like 
ceiling. The digits were based on an 1853 calculation by English mathematician William Shanks, which 



included an error beginning at the 528th digit. The error was detected 
in 1946 and corrected in 1949/ 208 ] 

In Carl Sagan's novel Contact it is suggested that the creator of the 
universe buried a message deep within the digits of Tl/ 209 ] The digits 
of 7T have also been incorporated into the lyrics of the song "Pi" from 
the album Aerial by Kate Bush/ 210 ] 

In the United States, Pi Day falls on 14 March (written 3/14 in the 
US style), and is popular among students/ 211 ] 7T and its digital 
representation are often used by self-described "math geeks" for 
inside jokes among mathematically and technologically minded 
groups. Several college cheers at the Massachusetts Institute of 
Technology include "3.14159"Z 212 ] Pi Day in 2015 was particularly 
significant because the date and time 3/14/15 9:26:53 reflected many 
more digits of pi/ 213 ] In parts of the world where dates are 
commonly noted in day/month/year format, 22 July represents "Pi 
Approximation Day," as 22/7 = 3.142857J 214 ^ 

During the 2011 auction for Nortel's portfolio of valuable technology patents, Google made a series of 
unusually specific bids based on mathematical and scientific constants, including 7i/ 215 ] 

In 1958 Albert Eagle proposed replacing 71 by T (tau), where T = tt/2, to simplify formulas/ 216 ] However, no 
other authors are known to use T in this way. Some people use a different value, T — 2/1 — 6. 28318. ,./ 217 ] 
arguing that T, as the number of radians in one turn, or as the ratio of a circle's circumference to its radius 
rather than its diameter, is more natural than 7T and simplifies many formulas/ 218 ® 19 ] Celebrations of this 
number, because it approximately equals 6.28, by making 28 June "Tau Day" and eating "twice the pie"/ 220 ] 
have been reported in the media. However, this use of I has not made its way into mainstream 
mathematics/ 221 ] 

In 1897, an amateur mathematician attempted to persuade the Indiana legislature to pass the Indiana Pi Bill, 
which described a method to square the circle and contained text that implied various incorrect values for 71, 
including 3.2. The bill is notorious as an attempt to establish a value of scientific constant by legislative fiat. 
The bill was passed by the Indiana House of Representatives, but rejected by the Senate, meaning it did not 
become a law/ 222 ] 

In computer culture 

In contemporary internet culture, individuals and organizations frequently pay homage to the number 7T. For 
instance, the computer scientist Donald Knuth let the version numbers of his program TeX approach 7T. The 
versions are 3, 3.1, 3.14, and so forth/ 223 ] 
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A pi pie. The circular shape of pie 
makes it a frequent subject of pi 
puns. 


Footnotes 


a. 


The precise integral that Weierstrass used was 7r = 



dx 

--. Remmert 2012, p. 148 

1 + x 2 K 


b. The polynomial shown is the first few terms of the Taylor series expansion of the sine function. 























c. Allegedly built so that the circle whose radius is equal to the height of the pyramid has a 
circumference equal to the perimeter of the base 
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